CONTROLLABILITY OF 3D LOW REYNOLDS SWIMMERS 
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Abstract. In this article, we consider a swimmer (i.e. a self-deformable body) immersed in a fluid, the flow 
of which is governed by the stationary Stokes equations. This model is relevant for studying the locomotion of 
microorganisms or micro robots for which the inertia eff'ects can be neglected. Our first main contribution is to prove 
that any such microswimmer has the ability to track, by performing a sequence of shape changes, any given trajectory 
in the fluid. We show that, in addition, this can be done by means of arbitrarily small body deformations that can be 
superimposed to any preassigned sequence of macro shape changes. Our second contribution is to prove that, when no 
macro deformations are prescribed, tracking is generically possible by means of shape changes obtained as a suitable 
combination of only four elementary deformations. Eventually, still considering finite dimensional deformations, we 
state results about the existence of optimal swimming strategies for a wide class of cost functionals. 
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1. Introduction. 

1.1. Context. Relevant models for the locomotion of microorganisms can be tracked back to 
the work of Taylor [TB] , Lighthill [TT| [TU] , and Childress [B] . Purcell explains in [H] that these 
sort of animals are the order of a micron in size and they move around with a typical speed of 30 
micron/sec. These data lead the flow regime to be characterized by a very small Reynolds number. 
For such swimmers, inertia effects play no role and the motion is entirely determined by the friction 
forces. 

In this article, the swimmer is modeled as a self deforming-body. By changing its shape, 
it set the surrounding fluid into motion and generates hydrodynamics forces used to propel and 
steer itself. We are interested in investigating whether the microswimmer is able to control its 
trajectory by means of appropriate shape deformations (as real microorganisms do). This question 
has already be tackled in some specific cases. Let us mention [T^ (the authors study the motion of 
infinite cylinders with various cross sections and the swimming of spheres undergoing infinitesimal 
shape variations) and [2 (in which the ID controllability of a swimmer made of three spheres is 
investigated). 

Our contribution to this question is several folds. First, we give a definitive answer to the control 
problem in the general case: the swimmer we consider has any shape at rest (obtained as the image 
by a diffeomorphism of the unit ball) and can undergo any kind of shape deformations (as long as 
they can also be obtained as images of the unit ball by diffeomorphisms) . With these settings, we 
prove that the dynamical system governing the swimmer's motion in the fluid is controllable in the 
following sense: for any prescribed trajectory (i.e. given positions and orientations of the swimmer at 
every moment) there exists a sequence of shape changes that make him swim arbitrarily close to this 
trajectory. A somewhat surprising additional result is that this can be done by means of arbitrarily 
small shape changes which can be superimposed to any preassigned macro deformations (this is 
called the ability of synchronized swimming in the sequel). Second, when no macro deformations 
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are prescribed (this is called freestyle swimming in the paper) , we prove that the ability of tracking 
any trajectory is possible by means of shape changes obtained as an appropriate combination of 
only four elementary deformations (satisfying some generic assumptions). Third, we state a result 
about the existence of optimal swimming. 

Notice that the paper follows the lines of 4, in which the authors study the controllability of a 
swimmer in a perfect fluid. 

1.2. Modeling. 

Kinematics. We assume that the swimmer is the only immersed body in the fluid and that 
the fluid-swimmer system fills the whole space, identified with R^. Two frames are required in the 
modeling: The first one € := (Ei, E2, E3) is fixed and Galilean and the second one e := (ei, 62, e^) 
is attached to the swimming body. At any moment, there exist a rotation matrix R e SO (3) and 
a vector r e R'^ such that, if X {Xi, X2, X^)* and x := {xi,X2,X3)* are the coordinates of a 
same vector in respectively £ and c, then the equality X — Rx + r holds. The matrix R is meant 
to give also the orientation of the swimmer. The rigid displacement of the swimmer, on a time 
interval [0,r] (T > 0), is thoroughly described by the functions t : [0,T] ^ R{t) e S0(3) and 
t : [0,T] 1—^ r(t) G R'^, which are the unknowns of our problem. Denoting their time derivatives 
by R and r, we can define the linear velocity v := (wi, W2, 1^3)* G R"^ and angular velocity vector 
r2 := (f2i, 5^2, ils)* e R'^ (both in e) by respectively v := R*r and O :~ R*R, where for every 
vector u := {ui,U2,U3)* G R^, u is the unique skew-symmetric matrix satisfying ux := u x x for 
every x € R'^. 

Shape Changes. Unless otherwise indicated, from now on all of the quantities will be ex- 
pressed in the body frame e. In our modeling, the domains occupied by the swimmer are images 
of the closed unit ball B by diffeomorphisms, isotopic the identity, and tending to the identity 
at infinity, i.e. having the form Id -I- ■(? where belongs to Dq(R"^) (the definitions of all of the 
function spaces are collected in the appendix, Section |A| . With these settings, the shape changes 
over a time interval [0, T] can be simply prescribed by means of functions t € [0, T] 1-^ -dt G Z?q(R'^) 
lying in M^^'^([0, T], Z3q(R'^)). Then, denoting 0t := Id + dt, the domain occupied by the swimmer 
at every time t > is the closed, bounded, connected set Bt '■— Ot{B) (keep in mind that we are 
working in the frame e) and wt ;= dtO{0~^) is the swimmer's Eulerian velocity of deformation. 
We shall denote S := dB the unit ball's boundary while St :— 0t(S) will stand for the body-fluid 
interface. The unit normal vector to St directed toward the interior of Bt is nt and the fluid fills 
the exterior open set Ft '■— R^ \ Bt- 

The Flovif. The flow is governed by the stationary Stokes equations. They read (in the body 
frame e): 

-AiAu -I- Vp = 0, V-u = inJ^t (t > 0), 

where /i is the viscosity, u the Eulerian velocity of the fluid and p the pressure. These equations 
have to be complemented with the no-slip boundary conditions: u = ft x x + v + wt on St. The 
linearity of these equations leads to introducing the elementary velocities and pressures (ui,pi) 
(i = 1,...,6) and {ud,pd)^ defined as the solutions to the Stokes equations with the boundary 
conditions = e^ x a; (i = 1,2,3), — ei_3 {i — 4,5,6) and — Wj on Ej. Then, the 
velocity u and the pressure p can be decomposed as u = J2i=i^i'^i + 12^=4^i-3'^i + and 
P = J2i=i^iPi +Si=4"i-3K + Pd- Notice that the pairs {ui,pi) {i = 1,...,6) and (urf,prf) are 
well-defined in the weighted Sobolev spaces (WQ^(J"f))^ x L'^{Tt) (see the Appendix, Section [€[). 
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Fig. 1.1. Kinematics of the model: The Galilean frame (£ := (Ej)i<j<3 and the body frame c := {e_,)i<j<3 
with Bj = REj (R a SO(3)j. Quantities are mostly expressed in the body frame. The domain of the body is Bt at 
the time t and Bt is the image of the unit ball B by a diffeomorphism 0f The open set Tt '■= R"^ \Bt is the domain 
of the fluid. The position of the swimmer is given by the vector r (in i£) and its orientation by R a SO (3). The 
vector V := _R*r is the translational velocity (in c). 



Dynamics. As already pointed out before, for microswimmers, the inertia effects are neglected 
in the modeling. Newton's laws reduce to J^^ T(u,p)nt x xda = (balance of angular momentum) 
and /j^^ T(u,p)nt da — (balance of linear momentum) where T{u,p) := 2fiD{u.) — pld is the stress 
tensor of the fluid, with D{u) := (Vu + Vu*)/2. The stress tensor is linear with respect to (u,p) so 
it can be decomposed into T(u,p) = X^Li ^i'^{ui,Pi) + Yf^=4Vi-3T{ui,pi) +T{ud,Pd)- In order to 
rewrite Newton's laws in a short compact form, we introduce the 6x6 matrix M(i) whose entries 
are 

^ f/j.^ e, • {T{uj,pj)nt X x)da = J^^{x x e,) • T{uj,pj)ntda (1 < i < 3, 1 < j < 6); 
\/s,e.-3-Tr(u,,p,)ntda (4 < z < 6, 1 < j < 6); 



and N(<), the vector of whose entries are 



/j.^ Gj • (T(urf,pd)nf X x)d(T = J^^{x X e^) • T{ud,Pd)ntda {1 < i < 3); 
/j. ei_3 • T(ud,Pd)ntdcr (4 < i < 6). 



With these settings, Newton's laws take the convenient form M(t)(r2, v)* + N(i) = 0. Upon an 
integration by parts, the entries of the matrix M{t) can be rewritten as Mij{s) := 2/i Jj-^ D{ui) : 
D{uj)dx, whence we deduce that M{t) is symmetric and positive definite. We infer that the 
swimming motion is governed by the equation: 

(^^) = -M(t)-iN(0, {0<t<T). (1.1a) 

To determine the rigid motion in the fixed frame 2;, Equation (1.1a) has to be supplemented with 
the ODE: 



dt[r ) ^ [rv 



iO<t<T), (1.1b) 



together with Cauchy data for i?(0) and r(0). At this point, we can identify the control as being 
the function t S [0,r] i— ?> G Dq(R^). Notice that the dependence of the dynamics in the control 
is strongly nonlinear. Indeed ??t describes the shape of the body and hence also the domain of the 
fluid in which are set the PDEs of the elementary velocity fields involved in the expressions of the 
matrices M{t) and N(t). 

Considering ( |1.1[ ), we deduce as a first nice result: 

Proposition 1.1. The dynamics of a micro swimmer is independent of the viscosity of the 
fluid. Or, in other words, the same shape changes produce the same rigid displacement, whatever 
the viscosity of the fluid is. 

Proof. Let {uj,pj) be an elementary solution (as defined in the modeling above) to the Stokes 
equations corresponding to the viscosity /i > 0, then {uj, {fl/^)pj) is the same elementary solution 
corresponding to the viscosity /i > 0. Since the Euler-Lagrange equation ( |1.1[ ) depends only on the 
Eulerian velocities Uj, the proof is completed. □ 

As a consequence of this Proposition wc will set /i = 1 in the sequel. 

Self-propelled constraints. For our model to be more realistic, the swimmer's shape changes, 
instead of being preassigned, should be resulting from the interactions between some internal forces 
and the hydrodynamical forces exerted by the fluid on the body's surface. To do so, the dynamics 



(1.1) should be supplemented with a set of equations (for instance PDEs of elasticity) allowing 
the shape changes to be computed from given internal forces. However, this would make the 
problem of locomotion much more involved and is beyond the scope of this paper. For weighted 
swimmers, this issue can be circumvented by adding constraints ensuring that the body's center of 
mass and moment of inertia are deformation invariant in the body frame. Unfortunately, massless 
microswimmers have no center of mass and their moment of inertia is always zero. 

To highlight the fact that constraints have still to be imposed to the shape changes for the 
control problem to make sense, consider the following result: 

Proposition 1.2. Let -d^ e W^'^{[0,T],D^(R^)) be two control functions such that := 
Id + d and 0^ := Id + i)^ differ up to a rigid displacement on the unit sphere (more precisely, 
for every t S [0,T], there exists {Q{t),s{t)) e S0(3) x such that (Q(0),s(0)) = (Id, 0) and 
= Q(i)0t|s + s{t)). Then, denoting by t G [0,T] ^ {R{t),r{t)) G S0(3) x a solution (if 
any) to System (1.1 1 with Cauchy data {Ro,ro) G S0(3) x R'^, we get that the function t G [0,T] t-^ 



{R^(t),r^t)) := {R{t)Q{t)*,r{t) - R{t)Q{t)*s{t)) G S0(3) x R^ is also a solution with the same 
Cauchy data but control . In particular R^ {t)0l {t) = R{t)0t+r{t) for all t G [0, T] (i.e. the 
swimmer's global motion is the same in both cases). 

Proof. If we denote by Ui(i) (i = 1, . . . , 6) (respectively uj(t)) the elementary velocity fields ob- 
tained with the control function ■& (respectively it can be verified that Ui{t, x) = Q(t)*ul{Q{t)x+ 
s{t)) for every t G [0, T] , every x € Tt and every i = 1, . . . , 6. We deduce that M(i) = Q(i)*Mt {t)Q{t) 
where the elements of M{t) (respectively M'^ (i)) have been computed with the elementary velocity 
fields Ui{t) (respectively u|(i)) and Q{t) G S0(6) is the bloc diagonal matrix dia.g{Q {t),Q{t)). On 
the other hand, denoting respectively by Wt(x) = dt0t{0t^{x)) and w|(x) = dt0\{0\''^ {x)) the 
boundary velocity of the swimmer in both cases, we get the relation: W((a;) + x(i) x a; + C,{t) = 
Q{t)*wl{Q{t)x + s{t)) for all t G [0,r], where x{t) ■= Q{t)*Q{t) and C,{t) := Q{t)*s{t). With 
obvious notation, we deduce that N(t) + M(t)(x(t), <(t))* = Q(t)*Nt(t). If we set now (V2,v)* := 
-M(i)N(t) and {Vl\^r^)* := -Mt(t)Nt(t), we get the identity (O^vt)* = Q(i)(f2 - x,v - C)*- 
It suffices to integrate this relation, taking into account that ((5(0), s(0)) = (Id, 0), to obtain the 
conclusion of the Proposition and to complete the proof. □ 
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If we apply this proposition with constant in time (the boundary of the swimmer is ©(S) 
at any time), we deduce that any shape change which reduces to a rigid deformation Q(t)x + s{t) 
on the swimmer's boundary (9(S) will produce a displacement {Q{t)* ,—Q{t)*s{t)). But if we 
compute the global motion of the swimmer, we obtain Q* {t){Q{t)0{x) + s{t)) — Q{t)*s{t) — 0{x) 
for every a; € E and every time t which means that the swimmer is actually motionless (the rigid 
deformation of the swimmer's boundary is counterbalanced by its rigid displacement). To prevent 
this from happening, we add the following constraints to the deformations (inspired by the so-called 
self-propelled constraints for weighted swimmers, see for instance [5]): 

[ Ot{x)da^0 (for aU t e [0,r]) and [ dt0t{x) x 0t{x) da = (for a.e. t e [0, T]). (1.2) 

About the existence of such deformations, we have in particular: 

Proposition 1.3. For every function d in VK^'^([0, T], L'o(R^)) such that f^0t=o{x)d(T 



0, there exists a function -d^ in W^'^{[0,T], DqCR!^)) satisfying (1.2) and an unique absolutely 
continuous rigid displacement t € [0,r] i— ((5(t),s(i)) € S0(3) x R-^ such that (5(0) = Id, s(0) — 
and eJis = {Q{t)0t + s{t))\^ for every t e [0,r]. 

In other words, the proposition tells us that any function of VF^'^([0, T], _Dg(R'^)) satisfying the 



first equality of ( |1.2[ ) at i = 0, can be made allowable (in the sense that it satisfies (1.2)) when 
composed with a suitable rigid displacement on the unit sphere. 

Proof. Define s{t) :— (l/47r) j.^0tda (an absolutely continuous function on [0,r]) and 0t := 
0t - s(t) for every t € [0,T]. The matrix J(i) := ||0t||^3ld - 0t ® 0td(J is always definite 
positive since (J(i)a;) ■ x — \\0t x xjl^adcr for all t S [0,T] and all x S R'^. We can then define 
x(t) := I{t)~^ j^dt0t X 0tda as a function of L^([0, T], R^). The absolutely continuous function 
t e [0,T] Q{t) e S0(3) is obtained by solving the ODE dtQ{t) = Q{t)x{t) with Cauchy data 
(5(0) — Id (we consider here a Caratheodory solution which is unique according to Gronwall's 
inequality). Then, we set s(t) := -Q{t)s{t) for aU t € [0,T]. The function 0t := Qit)0t + s{t) 



is in W^i^1([0,T],C1(R3)3), satisfies ([L2]) but does not take its values in Dl{R^) because 0t{x) = 
Q{t)x + s(t) + o(l) ^ X as |la:^||R3 — > +oo. Let and fi' be large balls such that UtG[o t] ^t{B) C 
and f2 C ri' and consider a cut-off function ^ valued in [0, 1] and such that ^ = 1 in SI and ^ = 
in R'^ . To complete the proof, define 0'^ as the flow associated with the Cauchy problem 
X{t, x) = £,{x)dtdt{x) + (1 - £,{x))dtMx). X(0, x) = 0t=o{x). □ 

Definition 1.4. We denote by A the non-empty closed subset o/ M^^'^([0, T], £)o(R^)) con- 



sisting of all of the functions verifying (1.2) 



1.3. Main results. The first result ensures the well posedness of System (1.1) and the conti- 
nuity of the input-output mapping: 

Proposition 1.5. For any T > 0, any function ?9 e W^'^{[0,T], DKR^)) (respectively of 
class , p = 1, . . . ,+oo,uj) and any initial data (i?(0),r(0)) £ S0(3) x R'^, System (1.1) admits 



a unique solution t e [0,T] H> (i?(t),r(i)) G S0(3) x R^ (in the sense of Caratheodory) absolutely 
continuous on [0,T] (respectively of class C^). 

Let {'&j)j>i C W^'^{[0,T], DlCR^)) be a sequence of controls converging to a function ■&. Let 
a pair (i?o,ro) G S0(3) x R^ be given and denote by t G [0,r] i-^ {R{t),f{t)) S S0(3) x R^ the 



solution in AC {[Q,T], SO {Z) y.'R^) to System (1.1) with control and Cauchy data {Rq^v a). Then. 



the unique solution {R^,r^) to System (1.1) with control 'd^ and Cauchy data (i?o,ro) converges in 
AC([0,T],SO(3) X R3) to (i?,f) asj-^+oo. 

We denote by M(3) the Banach space of the 3x3 matrices endowed with any matrix norm. 



The main result of this article addresses the controllability of System (1.1) 
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Theorem 1.6. (Synchronized Swimming) Assume that the following data are given: (i) A 
function d G A (the reference shape changes); (ii) A continuous function t e [0,r] H> {R{t),r(t)) G 
S0(3) X R"^ (the reference trajectory to he followed). Then, for any e > 0, there exists a function 
t G [0,T] t9t e Z?q(R'^) (the actual shape changes) in A, which can be chosen analytic, such that 

§0 = 4, ??T = ??T and sup,g[o^j,] (\\^t - Mcl,{m' + " ^WIIm(3) + II^W - r{t)\\n^] 



_ ^ __ ^ < £ 

where the function t G [0,T] i— >■ {R(t),r(t)) G S0(3) x R'^ is the unique solution to system (1.1) 



with initial data (i?(0),r(0)) = (-R(0),r(0)) and control 

This theorem teUs us that any 3D microswimmer undergoing approximately any prescribed 
shape changes can approximately track by swimming any given trajectory. It may seem surprising 
that the shape changes, which are supposed to be the control of our problem, can also be somehow 
preassigned. Actually, the trick is that they can only be approximately prescribed. We are going 
to show that arbitrarily small superimposed shape changes suffice for controlling the swimming 
motion. 

When no macro shape changes are preassigned we have: 

Theorem 1.7. (Freestyle Swimming) Assume that the following data are given: (i) A function 
d G Dq(R!^) such that J^Oda = (the reference shape at rest) (ii) A continuous function t G 
[0,r] ^ (R{t),f{t)) G S0(3) X R-^ (the reference trajectory). Then, for any e > there exists a 
function d G I?q(R^) (the actual shape at rest) such that (i) J^Oda = (ii) \\^ — ^I1_dJ(R3) < 
e and (Hi) for almost any A-uplet (Vi,...,V4) G (Cq(R'^)'^)'* satisfying J^Yidx — 0, J^0 x 
\ida — and /j^ x Vj dcr = (i,j — I,..., 4), there exists a function t G [0,T] i— >■ s{t) := 
{si{t), . . . , Si{t))* G R* (which can he chosen analytic) such that, using i9t + X]i=i Si{tyVi G 



Dq(R?) as control in the dynamics (l.ll, we ^et supjg jq — i?(t)||M(3) + ||r(<) — i"(t)||R3) < e 

where the function t G [0,T] i— >■ {R{t),v{t)) G S0(3) x R"^ is the unique solution to ODEs (1.1) with 



initial data (i?(0),r(0)) = (i?(0), f(0)). 

We claim in this Theorem that any 3D microswimmer (maybe up to an arbitrarily small mod- 
ification of its initial shape) is able to swim by means of allowable deformations (i.e. satisfying the 



constraints (1.2)) obtained as a suitable combination of almost any given four basic movements. 

If we still seek the control function i?i as a combination of a finite number of elementary 
deformations, i.e. in the form 

n 

^t=^ + J2s^m^, (1.3) 

where t e [0,r] h-> s{t) (si(t), . . . , s„(i))* G R" is in Li([0, T], R"), J^Odx = and (Vi,..., V„) G 
(Co(R^)'^)" is a fixed family of n vector fields satisfying J^Yidx — 0, J^0 x Yidx — and 
Vi X Vj dx = {i, j = 1, . . . ,n) we can state the following result: 
Theorem 1.8. (Existence of an optimal control) Let f : S0(3) xR^ x L)o(R^) x (Cq (R^)^) R 
be a continuous function, convex in the third variable and let K, he a compact ofH"^. Let (i?o, ^Oi'^o) 
and (i?i, ri, I?!) he two elements o/SO(3) x R-^ x DqCRP) such that there exists a control function "dt 



(i) having the form (1.3) with s(t) G K, for a.e. t G [0,T], (ii) satisfying 'dt=o = '&t=T = "^i and 



(Hi) steering the dynamics (1.1) from (i?o,ro) (at t — 0) to (i?i,ri) (at t — T). Then, among all 
of the control functions satisfying (i-iii), there exists an optimal control realizing the minimum 
of the cost 

T 

f(R(t),r{t),dt,dt^t)dt. 



The proofs of these results rely on the following leading ideas: First, we shall identify a set 
of parameters necessary to thoroughly characterize a swimmer and its way of swimming (these 



parameters are its shape and a finite number of basic movements, satisfying the constraints (1.2)). 
Any set of such parameters will be termed a swimmer signature (denoted SS in short). Then, the 
set of all of the SS will be shown to be an (infinite dimensional) analytic connected embedded 
submanifold of a Banach space. 

The second step of the reasoning will consist in proving that the swimmer's ability to track 
any given trajectory (while undergoing approximately any preassigned shape changes) is related to 
the vanishing of some analytic functions depending on the SS. These functions are connected to 
the determinant of some vector fields and their Lie brackets (we will invoke classical results of Geo- 
metric Control Theory). Eventually, by direct calculation, we will prove that at least one swimmer 
(corresponding to one particular SS) has this ability. An elementary property of analytic functions 
will eventually allow us to conclude that almost any SS (or equivalently any microswimmcr) has 
this property. 

Eventually, the existence of an optimal control in Theorcm |1.8| is a straightforward consequence 
of Filippov Theorem (see [U Chap. 10]) 

1.4. Outline of the paper. The next Section is dedicated to the notion of swimmer signature 
(definition and properties). In Section [s] we show that the matrix M(f) and the vector N(t) (in 
(1.1a)) are analytic functions in the SS (swimmer signature, seen as a variable) and in Section [4] 



we will restate the control problem in order to fit with the general framework of Geometric Control 
Theory. In the same Section, a particular case of swimmer will be shown to be controllable. In 
Section [5] the proof of the main results will be carried out. Section [6] contains some words of 
conclusion. Technical results and definitions are gathered in the appendix in order to make the 
paper more readable. 

2. Swimmer Signature. A swimmer signature is a set of parameters characterizing swimmers 
whose deformations consist in a combination of a finite number of basic movements. 

Definition 2.1. For any positive integer n, we denote C{n) the subset of DqCR"^) x {C^{RP)'^)"^ 
consisting of all of the pairs £:—{•&, V) such that, denoting :— Id + d and V '.— (Vi, . . . , V„), the 
following conditions hold (i) the set {V^js • e^, 1 < i < n, fc = 1, 2, 3} is a free family in C^(E) (ii) 
every pair (V, V) of elements of {O, Vi , . . . , V„} satisfies J^^ V da; = and Jj, V x V dx = 0. 

We call swimmer signature (SS in short) any element c of C{n). 

By definition, Dq{^^) is open in Cq(R'^)^ (see appendix, Section[A|. We deduce that for any 
c G C(n), the set {s := (si, . . . ,s„)* e R" : -i? + J27=i ^^'^i ^ ^o(R-^)} is open as weU in R" and 
we denote 5(c) its connected component containing s = 0. 

Definition 2.2. For any positive integer n, we call swimmer full signature (SFS in short) any 
pair c :— (c, s) such that c € C{n) and s £ S{c). We denote Ci?(n) the set of all of these pairs. 

Restatement of the problem in terms of swimmer signature (SS) and swimmer full 
signature (SFS). Pick a SS, c = (i?,V) G C{n) with V (Vi,...,V„) (for some integer n). 
Denote := Id + i? and for all s G S{c), 0^ := Id + z9 + X)"=i ■SiV^ (c := (c, s) G C{n) is hence a 
SFS). The body of the swimmer occupies the domain B := 0{B) at rest and Be ■— 0s{B) (for any 
s G S{c)) when swimming. Notice that within this construction, the shape changes on a time interval 
[0,r] (T > 0) are merely given through an absolutely continuous function t : [0,T] t-^ s{t) G S{c). 
If t G [0,T] I—)- s(t) G R" stands for its time derivative in L^([0, T], R"), the Lagrangian velocity at 
a point a: of S is X]r=i while the Eulerian velocity at a point x G /Be is X^iLi 
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with Wj(a;) :— 'Vi{0^^{x)). Due to assumption (ii) of Definition 2.1 the constraints (1.2) are 
automatically satisfied. 

The elementary fluid velocities and elementary pressure functions corresponding to the rigid 
motions depend only on the SFS. Therefore, they will be denoted in the sequel Uj(c) and Pi{c) to 
emphasize this dependence. The same remark holds for the matrix M(t) whose notation is turned 
into M[(c). The elementary velocity and pressure (udyPd) connected to the shape changes can be 
decomposed into = J2i=i^i^ii^) ^^'^ Pd = ^27=1 ^^''^ii^) respectively. In this sum, each pair 
(wi(c), 7ri(c)) solves the Stokes equations in Jc R-'^ \ with boundary conditions Wi(c) — 
on Ec :— dBc- 

Introducing the matrix N(c), whose elements are 



A^,(c) := 



J^Jx X e^) • T(wj(c),7rj(c))nd(T (1 < « < 3, 1 < j < n); 
/s e,_3 • T(wj (c),7rj(c))ndCT (1 < i < 6, 1 < j < n); 



(recall that the viscosity n can be chosen equal to 1), the dynamics (1.1a) can now be rewritten in 
the form: 



= -M(c)-iN(c)s, {0<t< T). (2.1) 



Let us focus on the properties of C(n) and Cpin). 

Theorem 2.3. For any positive integer n, the set C{n) is an analytic connected embedded 
submanifold of Cl{'R?f x (Co^(R3)3)" of codimension N := 3(n + 2){n + l)/2. 

The definition and the main properties of Banach space valued analytic functions are summa- 
rized in the article [TT^. 

Proof Foranyc:= (i?, V) G C(J(R^)^ x (C7(J(R^)^)", denote Vq := Id+i? and V := (Vi,...,V„). 
Then, define for = 0, 1, . . . , n, the functions : CHH^ f x (CpHR^)^)" ^ R3(»+i-fc) by Afe(c) := 

( /s It, ^ Vfc+i dx, ■ ■ • , /j] Vfe X V„ dx^ . Every function Afc is analytic and so is A := 

(Ao, . . . , A„)* : CoHR^)^ X (Coi(R3)3)» (N := 3(n + 2){n + l)/2). In order to prove that 

5c A (c) (the differential of A at the point c) is onto for any c E C{n), assume that there exist 
(n + 2){n + l)/2 vectors a{ G R^ (0 < i < j < n) such that: 



E 

i=0 



a, • {dcAic),c'') = 0, Vc" e C1(R3)3 X {C^{R^)f, (2.2) 



where a, := (a^, . . . , )* £ r3("+i-») (j ^ 0, ...,n) and c'* := (i9'',V'') G Co^(R^)^ x 
(Co^(R3)3)" with Vg := Id + ^f' and V'' := (V/, . . . , V^). Reorganizing the terms in ([2^, we 
obtain that: 



n „ k—1 n 

E / Vl!.[^a^^xV,+«^- ^ a^xv/dx^O. 



Since this identity has to be satisfied for any {-d^, V'') G Cq (R"^)-^ x (Cg (R'^))'^, we deduce that, for 
every fc = 0, . . . , n: 

fe— 1 n 

^a^^xV.ls + a^-^ ai. xV,|s = 0. (2.3) 
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Integrating this equality over E, we get that a'^ ~ {k ^ 0, . . . ,n). Taking into account Hypothesis 



(ii) of Definition 2.1 the identity (2.3 1 with fc = leads to chq — for every j = 1, . . . , n. There are 



no more terms involving Vq in the other equations and invoking again Hypothesis (ii) we eventually 
get = for 1 < I < j < n. So, equality ( |2.2| entails that a.i = for alH = 0, . . . , n and the 
mapping dcA{c) is indeed onto for all c G C{n). 

The linear space X — Ker9cA(c) is closed since A is analytic. Let Y be an algebraic supplement 
of X in Cg(R'^)'^ X (Co(R'^)'^)", and denote by Py the Hnear projection onto Y along X. A crucial 
observation is that the linear space Y is isomorphic to and hence it is finite dimensional and 
closed in Co1(R3)3 x {C^{R^)^)". Define the analytic mapping / : AT x F ^ R^ by fix,y) = 
A(c + X -\- y). The mapping dyf(0,0) = dcA{c) o Py being onto, the implicit function theorem 
(analytic version in Banach spaces, see [T7]) asserts that there exist an open neighborhood Oi of 
in X, an open neighborhood O2 of in Y, and an analytic mapping g : Oi Y such that 5(0) = 
and, for every {x,y) in Oi x 02, the two following assertions are equivalent: (i) f(x,y) = (or, in 
other words, c + x + y belongs to C(n)), and (ii) y = g{x). The analytic mapping g provides a local 
parameterization of C(n) in a neighborhood of c. 

In order to prove that C{n) is path-connected, consider two elements := {d^V^) and 
ct (#,V*) of C(n) and denote 0^ := Id + (V|,...,Vt) and 6>* := Id + i?*, 

(Vf,...,Vt). According to Definition A.2 ^'^(R^) is open and connected. This entails 



that it is always possible to find a continuous, piecewise linear path t : [0,1] ^ dt & £'q(R'^) 
such that 'dt=o = and 'dt=x — '9^ ■ We introduce = tg < ti < . . . < tfc = l, a subdivision 
of the interval [0,1] such that t z9t is linear on every subinterval [tj,tj+i] (j = 0, — 1) 

and we denote (9* := lA + dt, := dt=t,, & := Id + (j = 0, . . . , fc). Since Cq^R^)^ is an 
infinite dimensional Banach space, it is always possible to find by induction Wi, W2, . . . , W„ in 
C^iR^f such that (i) both families {Wi|s • e^, . . . , W„|s • e^, |e • e^, . . . , |s • e^, fc = 1, 2, 3} 
and {Wi|s • e^, . . . ,W„|e • e^, Vj|s • e^, . . . , V* |e • e^, fc = 1,2,3} are free in C^(R^) and (ii) 
for any pair of elements V, V, both picked in the same family, J^Ydx = 0, x Vdx = 

(for all j = 1, . . . , fc) and V x V'dx = 0. Define the function t e [0, 1] Vj e Cl{Il^f by 
Y\ (l-2i)V|+2tWi if < i < 1/2 and Vj := (2- 2i)W, + (2t- 1)V* if 1/2 < t < 1 and denote 
Vi :— (Vf , . . . , V") € (Cq (R'^)'^)". Eventually, a continuous function linking to is given by 
t e [0, l]^cte C{n) with ct (t?^ Vg^/a) if < i < 1/3, q := (i?3t-i, V1/2) if 1/3 < t < 2/3 and 
Ct-= (^9*,V3t/2_i/2) if 2/3<i< 1. □ 

We omit the proof of the following corollary, similar to that of the theorem above: 

Corollary 2.4. For any positive integer n, the setCp{n) is an analytic connected embedded 
submanifold ofC^(R^f x (C^^CR^f)-^ x R" of codimension N := 3(n + 2){n + l)/2. 

We denote by U the projection of C{n) onto D\{R^) defined by i7(c) = H for all c :— {1}, V) G 
C{n). The proof of the following corollary is a straightforward consequence of arguments already 
used in the proof of Theorem |2.3[ 

Corollary 2.5. For any positive integer n and for any 1) G n{C{n)), the section i7^^({i?}) 
is an embedded connected analytic submanifold of {d} x (Cg(R'^)'^)" (identified with (Cg (R'^)'^)"j 
of codimension 3n{n + 3)/2. 

3. Sensitivity Analysis of the Matrices M(c) and N(c). For any positive integers fc and 
I, we denote M(fc, I) the vector space of the matrices of size k x I (or simply M(fc) when / — fc). 

Theorem 3.1. For any positive integer n, the mappings c e Cpin) 1— > M(c) e M(6) and 
c e Cpin) ^ N(c) e M(6,n) are analytic. 

Let us begin with a preliminary lemma of which the statement requires introducing some 
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material. Thus, we denote F := Il^^\B (remember that B is the unit ball, E := dB and n is the unit 
normal to E directed toward the interior of B). For all ^ e £'q(R'^), we set S := Id + ^, := S{B), 
:= S{F) and Ej := H(E). We denote q := (C,>V), with W := {W\W^) C {C},{R^f)^, the 
elements of Q := Dl{R^) x (Co1(R3)3)2 and w| := W*(S-i) (i = 1,2). Finally, for every q e Q, 
we define: 

<?(q):= / D{ul^):D{ul)dx, (3.1) 

where, for every i = 1,2, there exists a function G L?{F^) such that the pair (Uq,Pq) € 
(^□^(J"^))^ X L'^{T^) solves the Stokes system: 

-Aujj + Vpjj = inJ-^, (3.2a) 
V-ujj = inJ-^, (3.2b) 
ujj = w^ onE^. (3.2c) 

The first equation has to be understood in the weak sense, namely: 

/ Vujj : Vvdx - / 33jj(V • v)dx = 0, Vv e {W^ {T^)f. (3.3) 

Recall that the function spaces are defined in the Appendix, Section \K\ 
Lemma 3.2. The mapping q G Q i-^ ^(l) € R- is analytic. 

Proof. We pull back equality (3.3) onto the domain F using the diffeomorphism S. We get: 

/ VUjjA^ : VVdx - [ P^M^ : VVda; = 0, VV e (Wp^ {F)f, (3.4a) 

where :^ < o 5 ", := pj^ o S', := det(VH), :== (yS*VS)-^J^ and := (VS"*)"! J^. 
Likewise, ( |3.2b|3.2c[ ) can be turned into: 

: VUjj = 0, in F, (3.4b) 

= W on E. (3.4c) 

We now claim that the mapping ^ e Dl{R^) h^- - Id e ^;g(R^M(3)) is analytic. Indeed, the 
mappings C e D^iR^) ^ VS*VS - Id e E°{R^,M{3)), A e E°{R^,M{S)) ^ (Id + A)-i - Id e 
£;g(R3,M(3)) and ^ G £'o(R^) >-> - 1 € Cg(R^) are analytic. Then, for i = 1,2, we define the 
analytic functions F' : Q x {Wt^{F)f x L^{F) {W^\F)f x L^{F) x (i?i/2(E))3 by: 

/(A5,U,.)-(B5,F, 
r'(q,U,P):= Be:VU 
V 7s(U-W^) 

where 7s : {W^{F))^ — > (//^/^(E))^ is the trace operator and 

(A^,U,V):= J VVA^-.VVdx, (V e {W\F) f , Y e {W^ {F)f 
{M^,P,Y):^ J PM^:VVdx, {P e L^F) , Y e {W^ {F) f) . 
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We wish now to apply the imphcit function theorem (analytic version in Banach spaces, as stated 
in [17) 1 to the analytic function F"^ . Observe however that we are only interested in the regularity 
result. Indeed, according to Proposition |C.2[ we already know that for all i = 1,2 and all q G Q, 
there exists a unique pair {U\,P\) € {W^{F)f x L'^{F) such that r*(q,U^,P^) = 0. For every 
q e Q, the partial derivative 9(c/^p)/^'(q, U^, P^) can be readily computed. Indeed, we have: 



'(A4,x, ■ 



7s(x) 



V(x,^)e(Woi(F)fxL2(s). (3.5) 



Let (f , ?7, g) be any element of {Wf^HF)f x (^p-^ x (iJ 1/2 (p))3^ 
(f,77, g), is equivalent to: 



The equation {d(u,p)r'ici, {U\, P^), (x, 



VxAj : VVdx- / ttBj : VVdx = (f,V) o , VV e (VF^ (i^))^, 

B5 : Vx = ?/, in F, 
X = g on S. 



According to Proposition C.2 there exists a unique solution (X;"") ^ (^o^l^))"^ ^ such that 

11x11 (WoHi^))^ + IklU"!^') ^ + ll'?IU"(F) + llg|l(Hi/2(E))3] where the constant Q > 

depends on ^ only. We infer that for every q G Q, d(ij^p)F^{(\, U^, P^) is a continuous isomorphism 
from {W^{F)f x L'^{F) onto (Wo"\F))3 x x {H^^^{^)f. The implicit function theorem 

applies and asserts that the mappings q G Q i-^- (Uq,Pq) G {Wq{F))^ x L'^{F) {i — 1,2) are 
analytic. 



To conclude the proof, it remains only to observe that the function <?(q) introduced in (3.1) 
can be rewritten, upon a change of variables as 



(q) = I |^(VC/iV^-i + (VC/^V^-i)*) : (VC/^V--^ + (Vt/^V--^) 



') dx. 



which is analytic as a composition of analytic functions. □ 
We can now give the proof of Theorem |3.1[ 

Proof. For any c := (c, s) G Cp{n), where c :— (i?,V), we apply the lemma with ^ := + 
1]"=! '^''^i ^'^d W\ G {ei X e^, i = 1,2,3} to get that the mapping c G Cpin) ^ M(c) G 
M(6) is analytic. To prove the analyticity of the elements of N(c), we apply the lemma again with 
e ^ + s^'V^, Wi G {e,; X S, e„ i = 1, 2, 3} and g {Vi, . . . , V„}. □ 

4. Control Problem. 

4.1. Controllable swimmer signature. Let us fix c G C{n) (for some positive integer n) 
and recall that S{c) is the connected open subspace of R" such that (c, s) G Cp{n). Introducing 
(fi, . . . ,f„) an ordered orthonormal basis of R", we can seek the function t G [0,T] s{t) G S{c) 
as the solution of the ODE s{t) = Y^"^i Ai(t)fi where the functions Ai : t G [0,r] 1— Xi{t) G R are 
the new controls, and rewrite once more the dynamics (2.1) as: 




-M(c,s)-iN(c,s)fj 
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(0 < t < T). 



(4.1) 



It is worth remarking that in this form, s is no more the control but a state variable and c G C{n) 



is a parameter of the dynamics. Considering (4.1), we are quite naturally led to introduce, for 



all c e CfH, the vector fields Xi(c) := -M(c)-iN(c)£, e R^ Y,(c) := (X|(c), X^(c), f,)* G 
TidS0(3) X R3 X R" (we have used here the notation X, := (X^, Xf )* e R^ x R^) and Zi{R, s) := 
7^i^Yi(c) e TrS0(3) x R^ x R" where TZr := diag(i?, i?. Id) G S0(6 + n) is a bloc diagonal matrix. 



The dynamics (4.1) and the ODE (1.1b) can be gathered into a unique differential system: 



= ^A,(t)Z^(i?,s), (0<t<T). (4.2) 



For every i = 1, . . . ,n, the function {R, r, s) e S0(3) x R^ x S{c) Zl{R, s) G TrS0(3) x R^ x R" 
can be seen as an analytic vector field (constant in r) on the analytic connected manifold M(c) := 
S0(3) X R^ X S{c). We denote C any element {R, r, s) £ M{c) and we define Z{c) as the family of 
vector fields (Z*)i<i<„ on M{c). 

Lemma 4.1. Let c be a SS fixed in C{n) (n a positive integer). If there exists C G M{c) such 
that dimLie^; -E(c) = 6 + n, then the orbit of Z(c) through any C G A4{c) is equal to the whole 
manifold M (c) . 

Proof. Rashevsky Chow Theorem (see T) applies: If Lie^Z(c) = T^A^(c) for all C G M.{c) 
(or more precisely, for all (i?, s) G S0(3) x S{c) since Z* does not depend on r) then the orbit of 
Z{c) through any point of M.{c) is equal to the whole manifold. Let us compute the Lie bracket 
[Z*(i?,s), Zi{R,s)] for l<i,j <n and {R,s) G S0(3) x S{c). We get: 



(X^ x X^)(c) 

[Zl{R, s), Zl{R, s)] = Ur I (XI X X2 ~ xj X X2)(c) | + Ui 



(a,.x|-5,^.xi)(c) 

(9,.x|-a,^.Xf)(c) I . (4.3) 




By induction, we can similarly prove that the Lie brackets of any order at any point ^ G M{c) 
have the same general form, namely the matrix TZr multiplied by an element of r(id,o,s)-'^(c). We 
deduce that the dimension of the Lie algebra at any point of M{c) depends only on s. According 
to the Orbit Theorem (see [1]), the dimension of the Lie algebra is constant along any orbit. But 
according to the particular form of the vector fields Z* (whose last n components form a basis of 
R"), the projection of any orbit on S{c) turns out to be the whole set S{c) (or, in other words, 
for any s G 5(c) and for any orbit, there is a point of the orbit for which the last component is s). 
Assume now that dimLie(;.Z(c) = 6 + n at some particular point :— {R* ,r*,s*) G M{c). Then, 
according to the Orbit Theorem, for any s G 5(c), there exists at least one point (iJ^, r^, s) G M{c) 
such that dimLie(i^^.rj,,s)-2^(c) = 6 + n. But since the dimension of the Lie algebra does not depend 
on the variables R and r, we conclude that dimLie(^Z{c) = 6 + n for all C G A^(c). □ 

Definition 4.2. We say that c, a SS in C{n) (for some integer n) is controllable if there exists 
C € A^(c) such that dimLieij-Z(c) = 6 + n. 

It is obvious that for a SS to be controllable, the integer n has to be larger or equal to 2. The 
following result is quite classical (a proof can be found in 5J): 

Proposition 4.3. Let c G C{n) (for some integer n) be controllable (with the usual notation 
c :— (i?, V), V :— (Vi, . . . , V„) and := + X]r=i -^i^i f'^''" s'very s G S{c) ). Then for any given 
continuous function t G [0,r] H' {R{t),Ylt),s{t)) G S0(3) x R'^ x S{c) and for any e > 0, there 
exist n functions Ai : [0, T] —> R (i = 1, . . . , such that: 



1. suptg[o,T] {\\Rit) - R{t)\\u{3) + \\r{t) - r(t)||R3 + ||i?j(t) - t?s(t) llci(R3)3) 
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2. R{T) = R{T), r(T) = r(T) and s{T) ^ s{T); 



where t G [0,T] i— > (_R(t), r(t), s(t)) e Al(c) is the unique solution to the ODE (4.2) with Cauchy 
data R{0) = R{0) e S0(3), r(0) = f(0) G R^, s(0) = s(0) G 5(c). 

Let us mention some other quite elementary properties that will be used later on: 

Proposition 4.4. 

1. If c := (i9,V) e C{n) (n > 2) is a controllable SS with V := (Vi,...,V„) G {C^{Il^)^)" 
then any c+ V+) G C(n + 1) such that V+ := (Vi, . . . , V„, V„+i) G {C^{Il^f)"+^ 
(for some V„+i G C™(R'^)'^J is a controllable SS as well. 

2. If c := (??, V) G C{n) (n > 2) is a controllable SS, then for any d* e {-d + s G 
iS(c)} the element c* := {■d*,V) belongs to C{n) and is a controllable SS as well. 

3. If c := V) G C{n) (n > 2) is a controllable SS, then all of the controllable SS in the 
section i7^^({i?}) form an open dense subset o/7T^^({i?}) (for the induced topology). 

4- If there exists a SS in C{n) for some n > 2 then, for any k > n, all of the controllable SS 
inC{k) form an open dense subset of C{k) (for the induced topology). 

Proof. The two first assertions are obvious so let us address directly the third point. Denote 
£fc {k positive integer) the set of all of the vectors fields on M.{c) obtained as Lie brackets of order 
lower or equal to k from elements of Z{c). Then, consider the determinants of all of the different 
families of 6 + jt. elements of £k as analytic functions in the variable V (the other variables and 
s = being fixed). Since c is controllable, there exist at least one k and one family of 6 + n elements 
in £k whose determinant is nonzero. According to Corollary |2.5| and basic properties of analytic 
functions (see [17 ), the determinant can vanish only in a closed subset with empty interior of the 
section iT^^({i9}) (for the induced topology). The proof of the last point is similar. □ 



4.2. Building a controllable swimmer signature. In this subsection, we are interested in 
computing the Lie brackets of first order [Z^(i?, s), Z^(i?, s)] at (i?, s) = (Id, 0), for a particular SS 
c := (Id, V) G C(4) (so the shape of the swimmer at rest is the unit ball). We make use of the usual 
notation V := (Vi, . . . , V4) (to be specified latter on), s = (si, . . . , S4) G S{c) and c :— (c, s). 

To carry out the aforementioned task, we introduce the classical spherical coordinates (p, a,/?) 
such that, for all x := (xi, X2, 2^3)* G R'^, a; 7^ 0, we have xi = pcos(a) sin(/3), X2 = £>sin(a) sin(/3) 
and Xz = gcos{(3). At each point (p, a, /3) we define the related local frame (eg,eQ,e^). For any 
n > 1, we call rigid spherical harmonics of degree — (n + 1) any function having the form: 



n 

(e,a,/3)^g-("+i) 7„y„,„(cos/3,a), (4.4) 

m— — n 



where j-n, ■ ■ ■ ,ln G R and Yn^m are the classical spherical harmonics of degree n G N and order 
m G {—n, . . . , n}. 

According to Lamb, [9] (one can also see the book of Happel and Brenner, [HI ch. 3.2, p. 
62]), the solution (u,p) of the Stokes equations around an immersed body of any shape can be 
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decomposed as follows (in the body frame): 



P^^P-in+i), (4.5b) 



+ 00 



where X-{n+i) and (/"-(n+i) are rigid spherical harmonics of degree — (n + 1). 

The functions p_(„+i), X_(„+i) and (/)_(„+i) (or more precisely the coefficients 7^ , k e {— n,...,n} 



arising in (4.4)) have to be determined in order to satisfy the boundary conditions on the surface 



on the body. This can be done following a method given in [Hj to which we refer for further details. 



The main interest of writing the solution of the Stokes equations in the form (4.5) is that the 
entries of the matrix M(c) and N(c) can be easily determined. 

Lemma 4.5. Let S be any smooth open bounded domain ofM.^ and denote J- :— R'^ \ S. Let 
(u,p) G (Wq{J-))^ X L'^{J-) be a solution to the Stokes equations given by (4.5) satisfying, for some 
no e N, x-(n+i) = 0-(n+i) = P-(n+i) = for all n > uq. 

For i — 1,...,6, let {ui,pi) G (VFq(J^))^ x L'^{J-) be the solution to the Stokes equations 
corresponding to the boundary condition u.i{x) = xxe.; ifi G {1, 2, 3} and\ii{x) — ei_3 ifi G {4, 5, 6} 
on dS . Then we have, 

KX r("'vl^»fe))' 

Proof. Let Ui be the rigid vector field defined by Ui(x) = a; x if « G {1, 2, 3} and VLi{x) — e^^a 
if i G {4,5,6}. Since u and are smooth, we have 2 /^Z?(u) : D{ui)dx — /g^T(u,p)ui . ndcr = 

T(u, p)ui -n dfj, where n is the normal to dS oriented towards the interior of S. Let i?(0, R) C 
be a ball centered at of radius R > such that S C -6(0, R) and denote ■= J^CiBlO, R). Using 
the Green formulae and the fact that for every i G {1, . . . , 6} we have D{u.i) (Vui + Vu*) /2 — 0, 
we obtain JggT{u,p)u.i ■ nda — — /g^^p T(u,p)u,; • ndcr, with n the normal to dJ-'n. oriented 



towards the exterior of Tji. Invoking the L^ orthogonality of the spherical harmonics, we get (4.6 1. 
□ 

When the body is specialized to be the unit sphere and the boundary conditions for u are x 



or ej {i,j = 1,2,3), the entries of the vectors in (4.6) are the elements of the matrix M(c, 0) and 



we get M(c, 0) = diag(87rld, 47rld). Similarly, if u = Vj (i = 1, . . . , 4) on the surface of the body, 
the entries of the vectors in ( |4.6[ ) turn out to be the elements of the matrix N(c, 0). Let now the 
vector fields be defined by Vi(g, a, /3) :— Vi{g, a, f3)eg for every i G {1, . . . , 4} with 

Vi{g,a,l3)=g-^''+^^n{Y3^i) (4.7a) 

V2{g,a,l3)^g-^''+'^^{Y3^i) (4.7b) 

V3{g,a,P)^g-^^+^^n{Y3^2) (4.7c) 

V4{g,a,P)^g-^^+'^n{Y4,2) (4.7d) 
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In this case, we get merely N(c, 0) = and hence Xi(c, 0) = {i ~ 1,...,4) in identity (4.3 1. 
Focusing now on the second term in the right hand side of (4.3), it remains to compute, for ah 
i,j — 1, . . . ,4 and c — (c, 0): 



9,.X,(c)-5„.X,(c) 



I(c)-i [(a,^.M(c)Xi(c) - a,^M(c)Xj(c)) 
I(c)-ik.N(c)f,-a,,N(c)f, 



(4.8) 



In particular, we need the expressions of the derivatives of the entries of the matrix N(c) with 
respect to s. 

Lemma 4.6. Let V e CqHR^)^ n C°°(R3)3 and wq G C°°(I])3 (recall that E is the boundary of 
the unit ball B and F ■.= R^\B). For every t small enough, we define Qt = Id + tV, Bt = Ot{B), 
:= dBt, Ft = -R^\Bt and w* - wo o 0^^ g C°°(I]t). 

Let also {ut,pt) and (uj,pj) G (VFj}(J-t))^ x L'^{Ft) (i — 1, . . . ,&) he the solutions to the Stokes 
problems in Ft with boundary conditions Ut = Wj, u\{x) = x x ei if i € {1,2,3} and uj(x) = e^-a 
if i (z {4,5,6} on Ej. Then we have 



d 
di 



[ D(uO : D{ul)dx) = [ Diu'o) : D{uU)dx, 



where Ug G {W^{F))'^ is solution of the homogeneous Stokes problem in F with the boundary 
condition 



-Vut=oV 



on S. 



(4.9) 



Proof. Since, for all t small, the solution Ut is smooth, according to [13 Theorem 4], the 
derivative of i i— > Uj at t = is solution of the homogeneous Stokes problem in F with boundary 
condition (4.9) (notice that the boundary condition is merely obtained by differentiating the equa- 
with respect to i at < = 0). Using the same argument as in the proof of Lemma 
D(uJ) da; = — Jgg(Q p) T(u(,pt)uj • ndtr. Differentiating with respect to 



tion Ut o 04 

45j we have 2/^^ D{ut) . ^^-t; - jaB{o,R) 
t and invoking the linearity of T and the Green formulae, we get the conclusion. □ 

1 



Applying Lemma 4.6 with the vector fields V.^ (i 
lengthy computations involving spherical harmonics: 



,4) defined in (4.7), we obtain after 
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One easily check now that dim (span {9s;N(c)fj — i9s^N(c)fi,l <i < j < 4}) — 6 and then 
dim (span {Z^(Id, 0), [Zj,(Id, 0), Z^(Id, 0)], 1 < A; < 4, 1 < i < j < 4}) =10 which is the dimension 
of S0(3) X R'' X S{c). It entails, according to the forth point of Proposition |4.4[ 

Proposition 4.7. For any integer n > A, the set of all the controllable SS is an open dense 
subset in C{n). 

5. Proofs of the Main Results. 



Proof of Proposition |l-5[ Let a control function -d be given in M^^'^([0, T], £'q(R^)) and 
denote Id + i3. With the notation of Lemma 3.2 at any time t the entries of the matrix M.{t) 
have the form <l>(q) with q := (i?,>V), W := {W^,W^), e {e^ x ©t, e„ i = 1,2,3} (j = 1,2). 
We deduce that t G [0, T] H' M{t) G M(3) is absolutely continuous. To get the expression of the 
elements of the vector N(i) we only have to modify which has to be equal to dt'dt- It entails 
that t e [0,T] N{t) e R'5 is in £^([0, T], R*"). Existence of solutions is now straightforward 
because t e [0,T ] ^ M(t)^iN(t) e R^ is in L^{[0,T],R^) and Caratheodory's existence theorem 



tion 



1.5 



denote by (il-' 
when the control is i?. As j 



applies to (1.1b I. Uniqueness derives from Gronwall's inequality. 

Let us address the stability result. With the same notation as in the statement of Proposi- 
v-')* the left hand side of identity (l.lal when the control is and (J7,v)* 

+00, it is clear that (n^v-?)* (l^,v)* in ^^([O, T], R^). Then, 

integrating (1.1b) between and t for any < i < T, we get the estimate \\R{t) — W {t)\\M{3) < 



/p ||i?(s) — i?-'(s)||M(3)|lf^(s)||R3 + ||Jl"'(s) — Jl(s)||R3ds. Applying Gronwall's inequality, we conclude 

that R in C([0, T], M(3)) as j — > +oo and we use again the ODE to prove that W R in 

^"'^([0, T], M(3)). Then, it is easy to obtain the convergence of r-' to r and to conclude the proof. 

Proof of Theorems |1.6| an d |1.7[ We shall focus on the proof of Theorem |1.6| because it 



will contain the proof of Theorem 



1.7 



For any integer n, we shall use the notation ||c| 



C(n) 



ll'^llciCR^)^ + Er=i l|V,|lci(R^)^ for alTc e Ci(R3)3 x (CoHR^)^)" with, as usual, c := V) and 
V:- (Vi,...,V„). 

Let e > and the functions t e [0, T] -dt & D^CR-^) and t e [0, T] i-> {R{t),r{t)) e S0(3) x R^ 
be given as in the statement of the theorem. According to Proposition |B.1[ we can assume that 
d £ C"^([0, T], _Dq(R'^)) n A because this space is a dense subspace of A. 

Step 1 (small initial jerking of the swimmer). In this step, we prove that the swimmer 
is able to modify slightly its shape in order to become controllable. Set := dt=o and V]^ := 
dt^t=o S C'o(R^)'^. According to the self-propelled constraints (1.2 1, it is always possible to find 
three elements V) {j = 2,3,4) in C|^(R^)^ such that the SS := (i?\ V^) belongs to C(4) (with 



(VI, 



(j 

,vl)). 



2,3,4) in Co1(R3)3 such that the SS 

guarantees that for any (5 > it is possible to find 



Then, Proposition 



a controllable SS in C(4), denoted by 



4.7 



i^\V^) where := {V\,. 



,vl) 



such that 



c ||c(4) < ((5 > is meant to be small an will be fixed later on). Moreover, 
can be chosen in such a way that there exists a smooth allowable function (i.e 
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we claim that 
satisfying 



pr2l )) t e [-1,0] ^ e_Dl{R^) such that = 1?^ and iJO^q = -d^ (i.e. the swimmer 

can modify its shape from into by self-deforming on time interval [—1,0]). Indeed, denote 
ci := 0\V^) e C(4) a controUable SS such that ||ci - c^\\c{4) be smah. Then define := 
+ (1 + t)(i9i - T^i) for every t £ [-1,0]. Since i:)ii(R3) is open, for - i9i||ci(R3)3 small 



enough, will remain in L)q(R'^) for all t E [—1, 0]. Then, Proposition 1.3 asserts that there exists 
a function Qo € ACQ-1, 0], S0(3)) and an allowable shape function ^ E W^'^i-l^O], Dl,(R^)) 
such that 6)? links 0^ (at t = -1) to some (at t = 0) satisfying 6>i|s = Qo(0)6>^|s- A 



careful reading of the proof of Proposition 1.3 allows noticing that ||Qo — Id||c([-i,o],M(3)) and 



¥t - ^'llM/i.Mhi,o],i3j(R3)) go to as II^Ji - t?i||ci(R3)3 goes to 0. Set now - Id, V} := 

Qo{0)^} {i = 1, . . . , 4) an d ob serve that the resulting SS satisfies the requirements. Furthermore, 
according to Proposition 1.5 [["^^ — '!^^||c^(r3)3 can always be made small enough for the control 



function 1?° to produce a rigid displacement t e [—1,0] i-> {Ro{t),ro{t)) G S0(3) x R'^ satisfying 
suptg[_i_o] iWMt) - ^(0)||m(3) + l|ro(i) - f(0)||R3 + - ^ici(R^)0 < Eventually, remark 
that this step of initial jerking performed on the time interval [toi^i] ■= [^IjO] can actually be 
carried out on a time interval arbitrarily short just by rescaling the time. 

Step 2 (building a continuous piecewise control function). Since the function dt^ is 
continuous on the compact set [0,r], it is uniformly continuous. For any v > 0, there exists 5^ > 
such that ||9t'i?t — dt''^t'\\c^(Ra)3 < v providing that |i — i'| < 6^. Then, we divide the time interval 
[0, T] into Q ^ ti < t2 < ■ ■ ■ < tk = T such that \tj+i - tj\ < for j = 1, . . . , fc - 1. For any 
t G [^1,^2], we have the estimate: 

11^, - (^1 + it- <l)V})||cl(R3)3 < - + {t- il)V})|lc.(R3)3 

+ W " ^'||C,1(R3)3 + (t - h)\\Vl - Vl|lcl(R3)3. 

On the one hand, we have, for all t g [^1,^2], ||i^t — {^^ + {t ^ ^i)Vi)||c1(r3)3 < I'lt — <i|. On 
the other hand, still for ti < t < t2 and if we assume that S^, < 1, we get — ^?^||i::1(R3)3 + 
(< — ti)||V} — V}||c;i(R3)3 < 6/2. We denote := '&t=t2- It is always possible to supplement 
:= dt-dt2 with vector fields (j = 2, . . . , 4) in such a way that := {¥, V^) be in C(4) with 
the obvious notation .= (v^, . . . , V|). We define i?^ ^1 + {t^ _ ti)Y\. For any ti<t < t2, 



Proposition 4.4 guarantees that the SS c\ := (i?^ + (^ — ^1)^1, V^) is controllable. In particular, 
for t = t2, there exists an integer k and a family of 10 vector field^ in £k (the set of all the 
Lie brackets of order lower or equal to k) such that the determinant of the family is nonzero. 
But this determinant can be thought of as an analytic function in V^. The set n~^{{'d^}) being 



an analytic connected submanifold of (Cg(R'^)^)^ (see Corollary 2.5), the determinant is nonzero 
everywhere on this set but maybe in a closed subset of empty interior (for the induced topology). 
Therefore, it is possible to find G {C^{R^)^y such that the SS := {d^, V^) is controllable and 
I|c'-c2||c(4) < {S/2 + uit2-h)) + S/4. 

By induction, we can build and (j — 1, 2, ... , k) such that (i) ||c-' — c^ || < '5/2 + X]i=2 + 
^{ti — ti^i) < S + vT and (ii) every c' is controllable. We choose 5 and v in such a way that 
6 + vT < e/4 and we define t : [0,T] 1— > -iJ^ G L'q(R'^) as continuous, piecewise afhne functions 
by -dt := -d^ + {t - tj)Y{ if t G [tj,tj+i] {j = l,...,fc - 1). Notice that for any t G [0,r], 
ll^t-^t|lc„MR3)3 <£/2. 



10 is the dimension of SO(3) X x S{cl^) 
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Definition 4.2 and Proposition 4.3 ensure that, on every interval [tj,tj+i] {j = 1, . . . , fc — 1), 



there exist four functions Xj : [tj,tj+i] R (i = 1, . . . , 4) such that the solution {Rj,rj,s^] 



[tj,tj+i] S0(3) X R3 X R4 to the ODE (4.2) with vector fields Z^, (i?j, s-?) and Cauchy data 



Riih) = i?o(0), ri(ii) = ro(0), R,{tj) = R{tj), v,{tj) = f(t,) {j = 2,...,k- 1) and s'itj) = 
(j = 1, . . . , fc — 1) satisfies: 

1- ^^VK[t,,t,^,] (\\m - RAt)\\MCi) + \\m - rj(t)||R3 + \\{}t - ^IWcHm^) < e/4 with := 

2. Rj{tj+i) = R{tj+i), rj{tj+i) = r{tj+i) and s^itj+i) = (i^+i - tj, 0,0,0)*. 
With these settings, the functions t e ^ e D},(R^), R : [-l,T] S0(3) and f : 

[-1, T] R3 defined by := ^l, R{t) := Rj{t) and f (t) := rj(t) if t e [ij, tj+i] (j = 0, . . . , fc - 1) 
are continuous, piecewise C^. 

Step 3 (smoothing the control function). We obtain a control function on [0, T] (still denoted 
by {)) by merely shifting/rescaling the time, from [— 1,T] onto [0,r]. Beforehand and as already 
mentioned, the first time interval [tQ,ti] [~1,0] could have been shortened as much as necessary 
for the estimate 



sup 

*e[o,T] 



to be true after the shifting/rescaling process. Then, we invoke Proposition B.l 



to conclude that there exists t e [0,r] n- i9f G Z3q(R'^) analytic, satisfying (1.2 1 and such that 



and Proposition 1.5 



sup - i?(t)||M(3) + \\r{t) - f(t)||R,3 + \\i}t- Mcam') < ^/2' 



*e[o,T] 



where (i?, r) : [0,T] ^ S0(3) x R^ is the solution to System ([TT]) with initial data (i?(0),r(0)) = 
(^(0),f(0)) and control The proof is then complete. 

6. Conclusion. In this paper, we have proved that every 3D microswimmer as the ability to 
swim (i.e. not only moving but tracking any given trajectory). Moreover, this can be achieved 
by means of arbitrarily small shape changes which can be superimposed to any preassigned macro 
deformation. When the shape changes are expressed as a finite combination of elementary defor- 
mations (and no macro shape changes are prescribed), we have shown that only four elementary 
deformations are needed for the swimmer to be able to track any trajectory. In this case and when 
the rate of shape changes (i.e. the velocity of deformations) is valued in a compact set, an optimal 
control exists for a wide variety of cost functionals. 

Appendix A. Function spaces. 

Classical function spaces. 

• For any open set ft C R'^ (included il — R'^), 2?(ri) is the space of the smooth (C°°) 
functions, compactly supported in Q. 

• For any open set C R'^ (included il = R"^), the set Cq{^1) is the completion of 'D{il) for 
the norm := sup^gf^ \u{x)\ + |j VM(a;)||R3. When il = R^, we get C^iR^) := {u G 
C^R^) : \u{x)\ and ||Vw(a;)||R3 ^ as ||a;||R3 ^ +cx)}. 

• The space Cq(R'^)"^ is the Banach space of all of the vector fields in R'^ whose every 
component belongs to Cq(R'^). 
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• For any Banach space E and any T > 0, C"([0, T], E) is the space of analytic functions on 
[0,r], valued in E. 

• Let now E be an open subset or an embedded submanifold of an Euclidean space and 
T > 0, then AC{[Q,T], E) consists in the absolutely continuous functions from [0,r] into 
E. It is endowed with the norm ||m|Uc([o,t],b) := supj^jp ^.j ||Mt||_E + Jq \\dtUt\\Edt. 

• CJJ^{ft, M{k)) (m an integer) is the Banach space of the functions of class C™ in valued 
in M(fc) (M(fc) stands for the Banach space of the k x k matrices, k a positive integer) and 
compactly supported in fl. 

• E™{n, M(fc)) stands for the connected component containing the zero function of the open 
subset {M e CS'{n,M{k)) : det(Id + M(x)) 7^ Va; G R^}. 

Lemma A.l. The set i)^(R^) {t? G CoHR^)^ s.t. Id + is a diffeomorphism of R^} is 
open in Cq (R^)^. 

Proof. The mapping e Cq^R^)^ ^ Si) := inf .^s^ (Id + V?9(x),e) • e e R (S^ stands for the 

unit 2 dimensional sphere) is well defined and continuous. For any i)q g Z3q(R'^), we have di}„ > 
and for all x,y ^ R'^ and e :— {y — x)/\y — x\ the following estimate holds: {y + ^{y) — x — d{x)) ■ e = 
\y — x\ /Q^(Id + Vi9(a; + te),e) • edi > \y — x\5i}. We deduce that Id + is one-to-one if d is close 
enough to -i^o- Further, still for close enough to i?o , Id -I- ^9 is a local diffeomorphism (according to 
the local inversion Theorem) and hence it is onto. □ 

Definition A. 2. We denote Dq(R'^) the connected component of Dq^R?) that contains the 
identically zero function. 

If ?? e Co(R'^)'^ is such that ||i?||ci;(R3)3 < li the local inversion Theorem and a fixed point 
argument ensure that Id -I- ?9 is a diffeomorphism so we deduce that Dq{R?) contains the unit 
ban of CoHR^)^. 

Sobolev spaces. 

• We define the weight function 9{x) := \J\ + |a;p (x G R^^) and the weighted Sobolev spaces: 

Wl{T) := {u e V\T) ■ Q-^u e L^{T)), (A.l) 
Wl (T) := {u e WliT) : 7s(h) = O}, (A.2) 

O 

where 7s : W^{T) 1— )■ is the classical trace operator. The dual space of {F) is 

• For any Banach space E, M^^'^([0, T], _E) is the Bochncr-Sobolev spaces (see for instance 
[T51 §7.1, page 187]) consisting in the functions w : [0,r] £' measurable and such that 
u and u' belong to L^([0,T],i?) (the derivative u' as to be understood in the sense of 
the distributions). It can be proved that W^''^{[0,T],E) is continuously embedded in 
C{[Q,T],E) and that u{t) = it(0) + M'(s)ds for all t G [0,T] and all u e W^-'^{[Q,T], E), 
where the integral is a Bochner integral (a generalization to Banach space valued func- 
tions of the Lebesgue integral). The space W^'^{[0,T],E) is endowed with the norm 
I|m||wi>i([o,t],b) := liu||c([o,T],£;) + /(f |lu'(s)|l£; ds. 

Appendix B. Control functions smoothing. 

Proposition B.l. For every e > Q and every d ^ A, there exists d ^ {[Q^T], Dl{'R?)) A 
such that \\§ ~ "i^ll wi'i([o,T],Di(R3)) < £ O'f^d ^t=a = i^t=o- In particular C"([0, T], Z?q(R^)) H A is 
dense in A. 
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Proof. Let i9 be in A. Since, C"([0, T],C^{R^)^) is dense in L^{[0, T],C^{R^f), we can always 
pick an element ( € C'^{[0,T],Cq{R?)^) which makes ||C~f^t'^llLi([o,T],C|;(R3)3) as small as required. 
Define for every t e [0,r] the analytic function dt = '!?t=o + /g C('5)'is, keeping in mind that the 
quantity — '&\\w'-''-{[o,t].dI(r^)) <^an be made arbitrarily small. Following the lines of the proof of 



PropositionlOj we define s{t) := (l/47r) 6>t da, Oj := Ot-s{t), the matrix J(t) := ||0||j2j^3ld- 
(9j(g)0jdcr (positive definite for every t e [0, T]) and x{t) :— J(t)^^ dtOjxOl da. Let us introduce 
aswellj(t) := Jj.\\0t\\^3id — Ot®0tda. Observing again that the quantity ||i9 — i9||y^^i.i([o,T].Di(R3)) 
can be arbitrarily small, we draw the same conclusion for ||s||vfi i([o,t],r3), then for ||J(i)~^ — 
J(i)"^llco([o,T],A/(3)) and finally for ||xIIli([o,t],r3)- Wc infer that HQ - Id||H/i,i([o,T],M(3)), the 
solution to the Cauchy problem dtQt = QtX{t) with initial data Qt=o ~ Id is arbitrarily small as 
weh. Then we set = Q{t)0\. At this point, 0* satisfies (Ol but (for t e [0,T]) is unlikely 



in DqCR.^) (because 0*{x) — )■ Q{t){x — s{t)) ^ x as ||a;||R3 — >■ +00). Notice however that for every 
smooth compactly supported function ^ : R'^ — R, the quantity ||C(t?*~'i?)||H'i i([o,T],_Di(R3)) can be 
made small. Let and Jl' be large balls such that lJte[o t] ^ti^) ^ ^ and ft C fl' and define ^ as a 
cut-off function valued in [0, 1] and such that ^ = 1 in and ^ = in R'^\f2'. To complete the proof, 
define as the fiow associated with the Cauchy problem X{t,x) — (,{x)dt'd'l{x) + {1 — (,{x))dt^t{x), 
X{0,x) = 0t=a{x). Indeed, - ^w^^{[Q.TlD},iR=^)) goes to as Uidti^t - 5t'?)llLi([o,T],ci(R3)3) 
goes to 0. □ 

Appendix C. Stokes Problem and Change of Variables. 

C.l. Well-posedness of the Stokes problem in an exterior domain. The following 
results that can be found in [7 : 

Theorem C.l. Let E be connected an Lipschitz continuous. Then, for any (f, h) G 
{Wq^{T))^ X L'^iT) X (iJi/2(I]))3, there exists a unique pair (u,p) G (W^iJ^))^ x L^(J") such 
that: 

-Au + Vp = f mT, (C.la) 
V • u = .9 m J", (C.lb) 
u = h on E. (C.lc) 

The solution has to be understood in the weak sense, namely: 



fvu:Vwdx-fpv-v)dx = (i,v) o , ywe(w^mf, 



(C.2a) 



V • u = g m J", (C.2b) 

7s(u) = h on E. (C.2c) 

Besides, there exists a constant Cjr > (depending on J- only) such that: 

\\^\\{w^{j^)Y + lblU"(.^) ^ ^A\\^\\(w-\fw + ll3lU"(.^) + l|h||(//i/2(s))3]. 



C.2. Change of variables. We denote, for all ^ G Dl{B,T{?), T := Id + ^, := det(Vr) 
and we define the matrices := (Vr*VT)-i and := (VT*)"! Jj. 
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Proposition C.2. //S is Lipschitz continuous, for all ^ e Dq^B^B.^) and for all (f, h) e 
{Wo\T)f X L^{T) X (iJi/2(jr))3 ^/jg following problem: 



VUeAeiVVdx- /P^B5:VVdx=(f,V) . 

B{ : VUj = g m J", 
7s(U^) = h on S, 



VVe«(^))3, (C.3a) 

(C.3b) 
(C.3c) 



admits a unique solution (Uj, P^) € (M^q^(J-'))'' x LP'[T). Moreover, there exists a constant C^{J-) > 
depending on J- and ^ only such that: 

\\~^i\\{w,]ir))^ + Wp^Wl^j^) < Q(-^)[l|f|!(i^-i(^))3 + hh^j^) + l|h||(i/i/2(s))3]- 



Proof. Let us introduce := Y(J-), := T(T,), :— g 
We denote by the distribution in (Wq"^(J^^))^^ defined by 



• T-^/{J^oT-^) and i^hoT- 



'(M/„-l(^))3x(Wol(^))3' 



This definition makes sense because there exist two constants ai(^) > (i = 1,2) such that 

ai{0\\v\\(wi(j-i)r ^ W ° ^\\{w-i{j^)r ^ oi2{0\Miwi{j-s)r ^ (W^o(-^?))^- Notice that 

when f is regular enough (i.e. can be identified with a function of {Ll^^(J^))'^) then we get merely 
:= f o T^^/{J^ o T^^). It is easy to check that, according to the properties of ^, the following 
mapping is a bicontinuous isomorphism: 

i?5 : {W^\T)f X L^{T) X {H^/^T)f ^ iW^\T^)f x L^{T^) x {H^/^{T^)f 

(f,.g,h)^(f^,g5,h^), 



Denote (u^,pj) = S^{f^, g^,h.^) the unique solution to the Stokes problem (C.2) in T^. The op- 
erator is hence a bicontinuous isomorphism mapping {Wq^{J^))^ x L^{J^) x {H^/^{J^))^ onto 
(M^g^(J^j))'^ X L^{J^(^). The following operator is a bicontinuous isomorphism as well: 

(v,g) ^ (V,Q)-(vor,9or). 



The solution to problem (C.3) is provided by the operator Tj :— o o and the following 
diagram commutes: 



(f,5,h) 



The proof is then completed. □ 
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